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Figs. 1 and 2. 1-Urge volcanic shield and genUe slopes covered with overlapping lava Hows; 2-centml domelike 
pan of some volcanic shields with many small volcanic hills and absence of expressed lava flows; 3-bright smooth 
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Fig. 3. Correlation of the mapped units shown in Figs. 1 and 2. 
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GEOID, TOPOGRAPHY, AND CONVECTION-DRIVEN 
CRUSTAL DEFORMATION ON VENUS. Mark Simons, 
Bradford H. Hager, and Sean C. Solomon, Department of Earth, 
Atmospheric, and Planetary Sciences, Massachusetts Institute of 
Technology, Cambridge MA 02139, USA. 

Introduction: High-resolution Magellan images and altim- 

etry of Venus reveal a wide range of styles and scales of surface 
deformation ( 1 ] that cannot readily be explained wi thin the classical 
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terrestrial plate tectonic paradigm. The high correlation of long- 
wavelength topography and gravity and the large apparent depths of 
compensation suggest that Venus lacks an upper-mantle low- 
viscosity zone [2-5]. A key difference between Earth and Venus 
may be the degree of coupling between the convecting mantle and 
the overlying lithosphere. Mantle flow should then have recogniz- 
able signatures in the relationships between surface topography, 
crustal deformation, and the observed gravity field [6,7], 

Model: We explore the effects of this coupling by means of a 

finite element modeling technique. The crust and mantle in these 
models are treated as viscous fluids. We solve both the equations of 
motion and the heat equation at evety time step using a modified 
version of the two-dimensional Cartesian finite-element program 
ConMan [8], A passive marker chain tracks the crust-mantle inter- 
face and permits variation in thecrustal buoyancy as well as specific 
crustal and mantle rheologies. These rheologies depend on compo- 
sition, temperature, and stress. In addition to the flow field, the 
stress field in the lithosphere, the surface topography, and the 
resulting geoid are readily calculated. The models presented here 
use an irregular finitc-elementmesh that is 50 elements high and 1 60 
elements wide. Our maximum resolution is in the 40-km-thick top 
layer, where each element is 2 km high and 5 km wide. In all, the 
mesh is 800 km in the horizontal dimension and 400 km in the 
vertical dimension. We impose free-sl ip boundary conditions on the 
top and side walls, with no flow through these walls. Flow at the 
bottom boundary is constrained to be vertical with no horizontal 
flow permitted. This last boundary condition gives a virtual 800 
km x 800 km box. The surface topography is calculated from the 
vert. cal stresses on the top wall of the box. Top and bottom 
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Fig. 1. Temperature contours at 130 m.y. for the (a) upwelling and 
(b) down welling models. Contour interval is 100°C; the top of the box is ai 
500°C. 

temperatures are fixed at 500°C and 1250°C respectively. Initially, 
we impose a linear temperature gradient across the lithosphere and 
set the rest of the mantle to be isothermal. We investigate two 
classes of models. Row in the first class is initiated with a sinusoidal 
temperature perturbation throughout the box. This results in mantle 
flow that is dominated by concentrated downwelling (Fig. lb). In 
the second class of models, flow is driven by a hot patch at the side 
of the box , thereby driving flow by concentrated upwelling (Fig . 1 a). 

Results: In all our models, convection produces horizontal 

compression al stresses in lithosphere above downwelling mantle 
and extensional stresses in lithosphere above upwelling mantle. As 
the convective vigor increases so does the magnitude of the stress. 
In models with constant-viscosity mantle overlain by a constant- 
viscosity crust, stress in the crust reaches values in excess of 100 
MPa in less than 100 m.y. We find that the rate of increase in 
compressive stress decreases with increasing crustal viscosity. This 
is because the stronger the crust, the more the development of the 
convective instability in the mantle driving the deformation is 
impeded. We also find that the magnitude of the peak compressive 
stress achieved above the downwelling increases with higher vis- 
cosities and/or with thinner initial crustal layers; the stronger the 
crustal lid, the more tractions from mantle convection are supported 
in the crust. Since force balance on the crust requires that shear 
traction integrated along the base be balanced by normal tractions 
integrated through its thickness, the thinner the crust, the larger the 
horizontal stresses. 

Both analytical models [9,10] and our numerical models of 
convection-induced crustal flow indicate that the amplitude and 
sign of the topography are highly time and rheology dependent. In 
genera], possible responses of the crust to mantle flow can be 


divided into three categories. The first involves little if any crustal 
flow, and topography results mainly from the transmission of 
normal tractions induced by density contrasts within the mantle. 
The second possible regime involves substantial crustal flow, with 
geologically rapid thickening over convective downwelling and 
thinning over convective upwelling. In this regime the effects of 
crustal thickness variations dominate the topography. A third pos- 
sible regime lies between the first two, with 4t in phase” deformation 
on short timescales and crustal flow on longer timescales. A strong 
mantle lithosphere tends to shield the crust from convective shear 
tractions, and topography results mainly from the transmission of 
convective normal tractions. A relatively weak lower crust facili- 
tates crustal deformation, and the isostatic effects of crustal thick- 
ness variations dominate the topography. 

Consideration of geoid to topography ratios (GTRs) can restrict 
which regime of crustal response is appropriate for Venus. The 
distribution of estimated GTRs for several highland regions on 
Venus is bimodal with two clusters around 10 and 25 m/km [5]. The 
positive correlation of long -wavelength gravity and topography 
implies that there are no major regions that have negative GTRs. To 
keep the GTRs positive, the geoid must follow the surface topogra- 
phy in sign at all times. Although the regime of negligible crustal 
flow can correctly predict both the sign and magnitude of the GTRs, 
it does not allow for the crustal deformation (i.e., flow) inferred 
from observations of tectonic features cm Venus. 

In contrast to the regime of negligible flow, that of time- 
dependent crustal flow generates topography that changes sign. 
Over a mantle downwelling, the topography is negative in the early 
stages of deformation and positive in the later stages of deformation; 
the converse holds over a mantle upwelling. During the transitional 
period the topography goes through zero, the geoid does not go 
through zero, and the GTR is unbounded; this singularity is not 
observed on Venus. An example of the evolution of the GTR for a 
model in the time-dependent regime is shown in Fig. 2. Note that 
before 90 m.y. both topography and geoid are negative, but at 90 
m.y. crustal thickening effects begin to dominate and the topogra- 
phy switches sign, causing the GTRs to first become unbounded and 
then negative. 

In the absence of a mechanism by which the sign of the geoid 
anomaly mimics that of the topography over a given upwelling or 
downwelling, only the regime of rapid crustal flow is plausible. In 
the case of mantle downwelling this would also require that the 



Fig. 2. Evolution of GTR wiihiime. At each time step, the GTR shown is the 
slope of the best fit line to the collection of geoid and topography values 
upward continued to 250 km above the surface nodes in the model. The result 
is for a downwelling model, with a crust of constant viscosity 50 times that of 
the mantle. In this model the si gn of the topography is time-dependent, but that 
of the geoid is not. The singularity al about 90 m.y. corresponds to the lime 
when topography crosses the zero datum. Note lhat ihe GTRs for this model 
are neither always positive nor always bounded. 
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lower mantle be more viscous than the upper mantle in order to 
produce the required positive geoid anomalies. This has already 
been shown to be true for the Earth, where the observed geoid highs 
over regions of mantle upwelling and regions of mantle down- 
welling are best explained by the presence of a strong lower mantle 
[1 1,12]. The large positive GTRs and the presence of large shield 
volcanos in certain highland regions on Venus, such as Beta Regio 
and Eistla Regio, are best explained as areas of mantle upwelling 
[5,13,14]. The regime of rapid cnistal flow predicts crustal thinning 
over the upwelling. However, the extensive partial melt and ensuing 
volcanism expected over such regions of mantle may outweigh the 
effects of crustal thinning on the surface topography and thus also 
yield postive GTRs [15]. 
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VENUS GRAVITY: SUMMARY AND COMING EVENTS. 
W. L. Sjogren, Jet Propulsion Laboratory, 4800 Oak Grove Drive, 
Pasadena CA 91109, USA. 

The first significant dataset to provide local measures of venu- 
sian gravity field variations was that acquired from the Pioneer 
Venus Orbiter(P VO) during the 1979-1981 period. These observa- 
tions were S-band Doppler radio signals from the orbiting space- 
craft received at Earth -based tracking stations. Early reductions of 
these data were performed using two quite different techniques. 
Estimates of the classical spherical harmonics were made to various 
degrees and orders up to 1 0 [ 1 ,2,3]. At that time, solutions of much 
higher degree and order were very difficult due to computer 
limitations. These reductions, because of low degree and order, 
revealed only the most prominent features with poor spatial resolu- 
tion and very reduced peak amplitudes. 

Another reduction technique was the line-of-sight acceleration 
mapping that had been used successfully for the Moon and Mars. 
This approach provided much more detail and revealed the high 
correlation of gravity with topography [4,5,6]. However, this tech- 
nique does not produce a global field as do the spherical harmonics. 
It provided a mapping of features from approximately 50°N to 25°S 
latitude for 360° of longitude. Other shortcomings were that the 
accelerations were at spacecraft altitude rather than at the surface 
and were not vertical accelerations; however, the reductions were 
quick and cheaply accomplished. Other efforts to analyze these data 
included local area reductions, where surface masses were esti- 
mated [7,8,9]. 

The computer revolution over the past 10 years has allowed new 
reductions with spherical harmonics. New fields up to degree and 
order fifty (2600 parameters) have been made [10,11,12], These 
fields now provide the best representation for any serious geophysi- 
cist doing quantitative modeling. There is now vertical gravity at the 
surface from a global model that carries all the requirements of 
dynamical consistency. There is one sizeable concern in that the 


resolution over the entire planet is not uniform. This is due to the 
Pioneer orbit, which had a high eccentricity, causing the high 
latitude regions of Venus to be poorly resolved. 

The Magellan (MGN) spacecraft, which went into orbit about 
Venus in August 1990, has returned Doppler data for gravity field 
reduction. However, because the high gain antenna was pointed at 
Venus for SAR mapping, no gravity data were acquired until the 
antenna was pointed back to Earth. This occurred at spacecraft 
altitudes higher than 2500 km, greatly reducing local gravity 
sensitivity. MGN has an eccentricity much smaller than PVO, so 
there is new information in the polar regions. Present reductions 
include two MGN circulations (486 days), which reduce uncertain- 
ties and produce somewhat better resolution. 

During March, April, and May 1992 new low-altitude data have 
been acquired from both PVO and MGN. PVO periapsis latitude has 
changed 27°, from 16°N to 11°S. These data will provide better 
definition in the southern hemisphere, particularly over Artemis. 
The MGN mission now acquires periapsis gravity data for one orbit 
out of eight (i.e., foregoes SAR mapping for one orbit/day). Since 
MGN has an X-band radio signal, the data quality is a factor of 10 
better than PVO. Only a small block of MGN data was acquired 
before its periapsis went into occultation May 16. Solar conjunction 
and periapsis occultation has also occurred for PVO. 

In September of 1992 MGN periapsis will exit occultation and 
its periapsis altitude will be lowered to approximately 170 km. 
Periapsis will be visible from Earth for a complete 360° longitude 
coverage period (243 days). This should be an excellent dataset, 
having low X-band data noise that in turn can be combined with the 
PVO dataset 

In December 1992 PVO will exit periapsis occultation and low- 
altitude data (~150 km) in the southern hemisphere will be acquired 
for about one month before PVO is lost due to the lack of fuel to 
maneuver to safe altitudes. 

In May 1993 there remains the possibility of aerobraking MGN 
into a circular orbit, thus allowing global uniform resolution gravity 
data to be acquired. One hopes that NASA has enough foresight to 
keep Magellan alive so this is a reality. It is anticipated that if this 
is done, harmonic solutions to degree and order 60-70 (5000 
parameters) will be produced. One could then compare similar 
features globally, resolve coronae and test many interior structure 
models. 
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DIFFERENT TYPES OF SMALL VOLCANOS ON VENUS. 
E. N. Slyuta 1 , I. V. Shalimov 2 , and A. M. Nikishin 2 , Vernadsky 
Institute, Russian Academy of Science, 117975 Moscow, Russia, 
2 Moscow University, Moscow, Russia. 

One of the studies of volcanic activity on Venus is the compari- 
son of that with the analogous volcanic activity on Earth. The 
preliminary report of such a comparison and description of a small 



